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PURPOSE  (U) 


The  Contractor  shall  furnish  the  necessary  personnel  and  facilities  for 
and,  in  accordance  with  instructions  issued  by  the  Scientific  Officer  or 
his  authorized  representative,  shall  conduct  an  applied  research  pro¬ 
gram  in  electronic  warfare  direction  finding  techniques: 

Against  HF  (1-30  MHz)  radio  signals  to  significantly  improve 
the  D/F  accuracy  achievable  with  today's  field  multipath  signals. 
Specifically,  to  devise  and  demonstrate  a  transportable  automated 
D/F  antenna  (and  associated  electronics  for  data  processing)  im¬ 
mune  to  interferring  signals  and  external  RF  noise.  In  particular, 
this  program  includes  the  following: 

1.  Study  and  conduct  socond-order  error  analysis  of  data  pro¬ 
cessing  methods  and  complete  the  antenna  concept  analysis. 

2.  Conduct  circuit  experiments  to  determine  the  feasibility  of 
goniometer,  broadband  receiver  and  broadband  data  processing 
concepts. 


ABSTRACT  (U) 


The  feasibility  of  utilization  of  a  short  dipole,  elevated  array  has  been 
demonstrated  on  ground-wave  signals  and  sky-wave  signals.  The 
feasibility  of  sky-wave  ranging  has  also  been  demonstrated.  Theoretical 
and  experimental  results  are  included. 

A  data  package  has  been  submitted  for  the  Interim  Replacement  for  the 
AN/TRD-21  Antenna  System. 
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HF  D/F  Concepts  for  Marine  Corps  AMPHIBIOUS  HF 
D/F  System 


Section  I 

CW  D/F  --  INTERIM  REPLACEMENT  FOR  AN/TRD-21  (U) 


Work  has  been  completed  on  the  elevated  H-Adcock  array  that  will 
replace  the  monopole  array  used  in  the  AN/TRD-21  D/F  System.  A 
data  package  was  submitted  to  the  Marine  Corps. 

1'.  D/F  BEARING  TESTS  (U) 

A  comparison  of  bearing  accuracy  for  the  AN/TRD-21  D/F  set  and 
the  interim  array  was  made  as  a  function  of  azimuth  angle.  This 
test  was  conducted  for  a  limited  number  of  bearing  cuts  and  more 
tests  should  be  conducted.  Table  1  shows  the  results  of  the  com¬ 
parison.  Figure  1  illustrates  the  number  of  occurrences  of  a 
bearing  angle  A  0 • 

2.  OCTANTAL  ERROR  (U) 

The  bearing  observed  on  the  azimuth  indicator  of  the  AN/TRD-21 


HF  D/F  System 

is 
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Russell,  R.  G. ,  et  al,  "Antenna  Study  for  D/F  Application",  DSD  R-192 
(Semi-Annual  Progress  Report),  1  July  1967  to  21  December  1967, 
Contract  NOnr  4798(00),  February  1968,  Section  I-B. 


TABLE  1  Comparison  of  Bearing  Accuracy  (U) 
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The  octantal  error  is  noted  to  be  a  function  of  element  spacing,  RF 
frequency,  and  true  azimuth  angle  of  the  RF  signal.  An  SDS  Digital 
Computer  was  used  to  plot  the  octantal  error  as  a  function  of  observed 
bearing  over  the  HF  range. 

The  results  of  the  plot  are  shown  in  Figure  2. 


Section  II 


SKY-WAVE  FEASIBILITY  STUDY  (U) 

INTRODUCTION  (U) 

Measuring  the  elevational  angle  of  a  sky-wave  to  determine  range 
information  has  been  shown  to  be  feasible  using  a  small  array  with 
an  automated  system.  The  range  of  operation  will  be  80  to  500  km 
with  a  range  accuracy  better  than  ±  10  per  cent  range. 

The  approach  used  was  to  measure  the  ionospheric  height  at  the 
receiver  location  and  to  measure  the  elevational  angle  of  arrival 
thus  leading  to  the  range  by  solving  the  triangle.  Over  the  range 
of  80  to  500  km,  the  flat- earth  approximation  is  valid.  Curved- 
earth  equations  were  used  on  all  experimental  ranging  data  pro¬ 
cessed  in  work  at  SURC2  .  The  problems  encountered  are: 

a.  Ionospheric  behavior 

b.  Measurement  requirements 

The  conclusions  reached  were: 

a.  Range  bracket  saves  problem 

b.  Additional  check  data  desirable 

The  contents  of  the  report  which  follows  will  cover: 

a.  Permissible  ionospheric  misbehavior 

b.  Ray  measurement  analysis  and  data 

c.  Planned  follow-on  work 


Boyko,  L. ,  Carmichael,  J. ,  CooperS.,  Cummins,  J. ,  Schuman, 
H. ,  "Antenna  Study  for  D/F  Application",  DSD  R- 193  ,  (Interim 
Report),  Contract  NOnr  4798(00),  October,  1968. 


IONOSPHERIC  BEHAVIOR  (U) 


a.  Introduction 

The  purpose  of  this  section,  is  to  present  and  evaluate  the  effects 
of  the  ionosphere  on  various  modes  of  sky-wave  propagation. 
Some  major  parameters  used  in  determining  and  describing 
the  characteristics  of  the  ionosphere  are: 

(1)  Layer  height 

(2)  Electron  density 

(3)  Layer  absorption 

(4)  Time  variations 

(a)  Diurnal 

(b)  Seasonal 

(5)  Gross  spatial  variations 

(a)  Polar 

(b)  Temperature 

(c)  Equatorial 

All  factors  are  markedly  influenced  by  solar  activity. 

Figure  4  shows  the  ionospheric  layers  of  most  interest  con¬ 
cerning  ionospheric  propagation.  A  flat-earth  description 
usually  is  adequate  at  short  ranges. 

Figure  5  displays  some  of  the  notation  used  throughout  this 
section.  Unless  otherwise  specified  the  letter  "F"  will  refer 
to  the  F3  layer. 

In  this  report  reliance  will  be  placed  on  the  ITSA  prediction 
technique3  which  primarily  considers  HF  reflection  from  the 
E  and  F  layers  and  absorption  from  the  D  layer. 

3  Lucas,  D.  L.  ,  Haydon,  G.  W. ,  and  Staff,  "Predicting  Statistical 
Performance  Indexes  for  High  Frequency  Ionospheric  Telecom¬ 
munication  Systems",  Environmental  Sciences  Services  and 
Aeronomy,  Frequency  Utilization  Section,  ESSA  Technical  Report, 
August,  1966. 
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FIGURE  3  Short  Base -Line  Operation  (U) 


A  55  66 


FIGURE  5  Geometry  and  Notation  Showing  One -Hop 
Propagation  (U) 
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b.  Tools  Available  to  Operator 

The  following  information  is  available  to  the  operator: 

(1)  Knowledge  of  latitude 

(a)  Polar 

(b)  Equatorial 

(c)  Temperature 

(2)  Knowledge  of  temporal  information 

(a)  Seasonal 

(b)  Diurnal 

(c)  Sunspot  cycle 

(3)  Real  time  information 

(a)  Recognition  of  spurious  effects 

(b)  Rejection  versus5  compensation  and  use 

c.  Ionospheric  Error  Analysis 

For  short  ranges  (R  <  50  km),  a  flat- earth  approximation  is 
acceptable  for  most  applications.  However,  as  a  check  the 
greater  part  of  the  analysis  summarized  in  this  section  was 
carried  out  for  a  spherical  earth.  Comparison  with  results 
for  a  flat  earth  indicated  that  use  of  the  approximation  was 
valid. 

In  Figure  6  the  variation  of  zenith  angle  (0)  with  range  (R)  is 
shown  for  two  virtual  heights  and  two  modes,  one  and  two  hop 
per  layer.  Note  that  the  zenith  angle  as  used  here  is  (90°  -  p) 
where  0  is  the  launch  angle.  To  represent  the  F2  layer,  a 
virtual  height  of  300  km  is  assumed  and  for  the  E  layer,  100 
km.  This  graph  gives  the  variation  of  zenith  angle  with  distance 
for  ranges  in  the  area  of  8  £  R  5  50  km. 

Since  current  information  on  the  state  of  the  ionosphere  is 
essential,  the  final  system  considers  the  availability  of  an  on-line 
ionosonde.  This  is  necessary  in  order  to  determine  the  virtual 
height  corresponding  to  the  equivalent  vertical  incidence  frequency 
(f^).  (f^  is  found  from  the  operating  frequency  by  the  "secant 

i  law"  where  i  is  the  angle  of  incidence  of  the  ray  upon  the  layer. ) 
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For  planning  purposes  characteristics  at  the  path  reflection 
point  are  frequently  obtained  from  the  ionospheric  predictions. 
While  these  serve  a  useful  .purpose,  it  must  be  recalled  that 
they  represent  the  general  climatology  of  the  ionosphere  for 
the  location  and  time  considered.  It  is  far  superior  to  employ 
current  knowledge  of  the  state  of  the  ionosphere  at  the  path 
midpoint.  In  a  practical  situation,  these  observations  may  be 
difficult  to  attain.  However,  it  usually  is  possible  to  make 
vertical  incidence  observation  on  a  noninterfering  basis  at  the 
station  site.  The  problem  then  becomes  one  of  extrapolating 
the  data  from  overhead  to  the  reflection  point.  On  a  general 
basis  the  greater  the  distance  the  poorer  the  extrapolation. 
Since  the  E  and  Fx  layers  so  closely  follow  the  behavior  in¬ 
dicated  by  the  ideal  Chapman  layer4  ,  their  characteristics  may 
be  projected  with  fair  accuracy  to  considerable  distances.  Un¬ 
fortunately  the  F2  layer  is  not  so  well  behaved  and  it  becomes 
difficult  to  define  an  accurate  extrapolation  function.  Practice 
has  shown,  however,  that  the  vertical  incidence  soundings 
made  at  one  point  hold  fairly  well  to  distances  of  about  150  km 
without  too  much  error  while  the  error  will  increase  with 
increasing  distance  from  the  measuring  point. 

Note  that  the  present  study  is  concerned  with  ranges  of  the 
order  of  80  £  R  £  500  km  implying  to  a  first  approximation 
that  reflection  points  occur  between  40  to  250  km  from  the 
site.  For  these  distances  information  supplied  by  an  ionosonde 
located  at  the  site  should  not  differ  too  much  from  that  at  the 
reflection  point  and  should  therefore  offer  considerable  im¬ 
provement  over  the  predictions.  For  these  distances  the 
sounder  data  may  be  employed  directly,  or  simple  correc¬ 
tions  for  distance  may  be  introduced. 

Figure  7  shows  the  effect  of  ionospheric  tilts  in  producing 
errors  in  range  (6R).  Let  a  =  tilt  of  the  layer  downward  tc 
the  receiver  with  respect  to  the  great  circle  between  the  two 
sites.  From  Figure  7 

03  =  9X  +  2a  for  ionospheric  tilt  case. 

4  Davies,  Kenneth,  "Ionospheric  Radio  Propagation",  National 
Bureau  of  Standards,  pg,  18. 


18 


FIGURE  7  Distortion  of  Ray  Path  Due  to  Tilt  (U) 


Then 


5R  =  h(tan  (91  +  2  a)  -  tan01  )  (flat  earth) 

The  range,  R,  varies  as  a  function  of  virtual  height,  h;  zenith 
angle,  0;  and  tilt  angle,  a,  as  follows: 


R 

R 

where 

a0 

t 


h  tan  0  +  h  tan  (1  +2 a)  (flat  earth) 


-l 


a0  rr  +  a0  0  -  2a0  a  -  a0  sin  j  sin  (0  +  2a) 


.  sin  (0  -  sin  t)  sin  2a  ’!  .  -i  .  , 

+  - - -  -  2a  sin  t  (curved 

earth) 


radius  of  earth  =  (6370  km) 


a  +  h 


sin© 


The  curved- earth  equation  was  used  in  constructing  the  curves 
of  Figure  8.  In  Figure  8  resulting  errors  in  range 
6  R 

|  —  |  x  100%  due  to  various  tilt  angles  (a)  are  shown.  Com- 

putations  were  carried  out  with  spherical  earth  approximations, 
and  representative  heights  of  100  and  300  kms  were  used  as 
additional  parameters.  Although  curves  for  only  positive  tilt 
angles  (a  =  0.  5°  ,  1°  ,  2°  )  are  plotted,  negative  tilt  angles 
were  found  to  have  almost  identical  corresponding  curves; 
e.  g.  ,  the  range  error  for  these  short  distances  due  to  a  =  -1° 
is  nearly  the  same  as  that  corresponding  to  a  =  1°  . 

Under  the  criterion  of  a  tolerable  error  in  range  of  10  per 
cent,  only  the  combined  conditions  of  short  range  and  high 
heights  appear  worthy  of  our  concern. 

Bramley  and  Ross6  have  shown  that  typical  tilt  (rms)  for  the 
F  layer  vary  from  ±  1°  to'±  2°  and  are  even  less  for  the  E 
layer.  Their  durations  may  vary  from  minutes  to  hours  with 

6 

Bramley,  E.  N. ,  Ross,  W. ,  "Measurement  of  the  Direction 
of  Arrival  of  Short  Radio  Waves  Reflected  at  the  Ionosphere", 
Proceedings  of  the  Royal  Society,  A,  220,  1951. 
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FIGURE  8  Per  Cent  Error  in  Range  Due  to  Tilt  versus 
Range  for  Various  Tilt  Angles  (U) 


perhaps  a  mean  of  5  to  30  minutes.  They  arise  from  the 
undulating  structure  and  movements  of  the  layers. 

Another  consideration  is  that  of  Sporadic  E(Eg).  Eg  is  a 
heavily  ionized  layer  near  the  normal  E  layer  in  height.  It 
appears  in  scattered,  isolated  instances  which  makes  pre¬ 
diction,  with  any  degree  of  accuracy,  considerably  difficult. 

Reflection  from  Sporadic  E  may  "blanket"  or  prevent  reflection 
from  the  normal  F2  layer.  This  blanketing  effect  extends  up 
to  the  equivalent  vertical  incidence  blanketing  frequency  of  Eg. 
Under  conditions  where  blanketing  Ec  occurs,  azimuthal  accu- 
racies  for  one-hop  propagation  may  exceed  those  obtained 
from  any  other  mode.  This  statement,  however,  neglects  pos¬ 
sible  tilts  of  the  Es  layer6*7  and  the  possibility  of  edge  or  side 
reflection  from  the  Es  cloud.  Also,  the  presence  of  Es  may 
introduce  other  complications.  If  a  localized  clouds  of  limited 
extent  is  present,  complex  modes  (e.  g.  ,  M,  N,  etc. )  may 
appear  making  a  determination  of  range  from  a  knowledge  of 
the  vertical  arrival  angle  difficult.  Note  that  when  the 
ranges  are  relatively  small  ameliorating  factors  appear.  For 
ranges  less  than  about  500  km,  the  probability  is  fairly  high 
that  only  single  or  double-hop  reflection  can  occur.  Measured 
launch  angles  when  double  hop  is  present  always  exceed  those 
for  single  hop  (see  Figure  6). 

It  is  possible  to  recognize  the  presence  of  double  hop  FJg 
propagation  fairly  well  at  least  during  the  summer  months. 

At  this  time,  in  the  northern  hemisphere 

(1)  values  of  f0  F2  are  relatively  low  and 

(2)  values  of  f0  Es  are  relatively  high. 

Also,  Sporadic  E  occurs  with  greatest  frequency  during  the 
period  0800-2000  LST.  If  during  this  period  high  vertical  angle 

g. - 

Harris,  R.  D. ,  Jensen,  D.  R. ,  Clark,  C. ,  "Analysis  of  Back- 
scatter  Studies",  Final  Report,  Utah  State  University,  Logan 
1965. 

7 

Gerson,  N.  C. Geddes,  W.H. ,  "Propagation  Data  Analysis", 
Second  Quarterly  Progress  Report,  1  October  1966  through 
30  December  1966,  Contract  DA18-  119-AMC-0 1952(X),  Task 
No.  14,  Electronics  Research  Laboratory,  Syracuse  University 
Research  Corporation,  Syracuse,  New  York,  27  February  1967. 
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reception  is  noted,  it  would  appear  probably  that  double  hop 
Es  propagation  had  occurred.  While  this  assumption  is  prob¬ 
ably  valid  in  the  great  percentage  of  cases,  information 
from  an  on-line  vertical  ionosonde  would  greatly  assist  in 
defining  the  modes  involved. 

Also,  due  to  these  characteristics  of  Es  during  the  summer 
months  and  because  the  ranges  under  consideration  are 
relatively  short  (<  500  km),  a  receiving  operator  can  be 
fairly  well  assured  that  if  an  oblique  signal  is  received  and  an 
Es  layer  is  determined  to  be  directly  overhead  then  the  in¬ 
coming  wave  is  traveling  only  via  the  Eg  layer  and  not  by 
either  the  Fs  layer  or  a  combination  of  both  F8  and  Es 
(N-type  mode).  This  conclusion  is  arrived  at  by  recognizing 
that  the  short  ranges  of  operation  require  for  F3  propagation 
large  elevation  angles  (or  small  zenith  angles  (Figure  7)),  and 
so  the  low  Fs  critical  frequency  approximates  the  MUF  of  the 
F3  layer.  Hence,  as  long  as  the  operating  frequency  is  not 
too  low,  an  absence  of  Fg  layer  propagation  can  usually  be 
assumed.  This  requires  any  propagation  to  be  via  the  Es 
layer. 

A  quick  check  can  be  performed  to  determine  with  good  prob¬ 
ability  from  which  layer  (Eg  or  F2  )  the  signal  is  reflecting. 

If  the  operating  frequency  f  satisfies 

f  <  4  E„  sec  i 
o  s 

where  f0  Eg  is  the  vertical  incidence  critical  frequency  of  the 
Sporadic  E  layer  and  is  is  the  apparent  angle  of  incidence  that 
the  ray  makes  at  this  layer  (if  it.  reflects  from  it),  then  the 
Es  and  not  the  F2  is  most  likely  the  layer  used  for  propagation. 

The  angle  ig  can  be  approximated  from  geometric  considerations 
by  the  relation 

cos  3 

h'  +  a0 
s  0 

radius  of  earth 

virtual  height  of  Es  (may  be  approximated  to  100  km) 
angle  of  arrival 


sin  i 


s 


where 

*o  = 

h'  = 
s 

0  = 
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An  on-line  computer  may  be  available  to  facilitate  the  com¬ 
putation. 

In  Situ  Tilt  Check 

Checking  for  tilt  conditions  can  be  made  by  making  known 
station  comparisons  for  a  reliability  assessment.  Another 
technique  is  to  look  for  time  variations  in  range  (especially 
day/night)  for  reliability  assessment  and  correction. 

Ionosonde  Prediction  Studies  Procedure 

A  comprehensive  study  has  recently  been  completed  demonstrating 
which  modes  of  operation  can  be  expected  and  the  relative  power 
loss  between  them.  Only  the  common  modes  were  considered: 
one  and  two  hop  E  and  F2  reflections  (see  Figure  9).  The 
parameters  considered  were  diurnal  variation  (hourly), 
seasonal  variation  (represented  by  the  months  December, 

March,  June  and  September),  geographical  location  (arctic, 
temperate,  and  tropic),  operating  frequencies  (2-12  MHz), 
and  transmitter  to  receiver  ranges  of  50,  250,  and  500  km. 

The  year  1967  was  chosen  as  a  model  year  for  acquiring  data  on 
which  to  base  the  computations.  -  This  is  primarily  because 
this  was  the  most  recent  year  of  which  an  entire  set  of  ITSA 
fQ  F2  and  M(3000)  F2  factor  prediction  data  cards  were  avail¬ 
able.  The  average  sunspot  number  (12-month  moving  average 
as  observed  at  Zurich,  Switzerland)  for  the  year  1967  was  96. 

The  procedures  used  for  determining  the  occurrence  of  various 
modes  of  propagation  are  basically  variations  of  the  procedures 
outlined  in  ITSA-13. 

In  the  employ  of  the  ITSA-1  prediction  technique,  the  following 
assumptions  were  utilized: 

(1)  Lower  level  refraction  was  ignored. 

(2)  Shimazoki's  formula8  for  computing  the  maximum  ioniza¬ 
tion  height  of  the  F2  layer  was  assumed  adequate. 

8 

Shimazoki,  F.  ,  "World-Wide  Daily  Variations  in  the  Height 
of  the  Maximum  Electron  Density  of  the  Ionospheric  F2 
Layer",  J.  Radio  Res.  Labs,  2,  No.  7,  Japan,  January 
1955. 


(3)’  Parabolic  layer  theory  was  used  throughout. 


f.  Results 

Figures  10  through  15  show  the  results  of  the  predictive  model 
study.  It  is  noted  that  usually  only  two  modes  are  effective 
within  500  km.  Equally  important  is  the  fact  that  th§  stronger 
mode,  usually  the  IF,  exceeds  the  others  by  10  dB  in  most  cases 
(see  Table  2).  This  latter  fact  will  be  used  to  advantage  in  the 
signal  processing  which  is  described  herein. 

g.  Conclusions 

The  results  of  this  study  indicate  the  following  conclusions 
should  be  made: 

(1)  Ionosphere  well-behaved 

(2)  Anomalies  detectable,  sometimes  correctable 

(3)  Normally  two  rays  present 

(4)  One  ray  generally  stronger  by  10  dB  to  30  dB 

(5)  For  80  km  £  R  £  500  km,  zenith  angles  ( 9 )  are  <  45 
degrees 


3.  WAVE  INTERFERENCE  (.U) 

The  following  discussion  uses  the  system  of  coordinates  shown  in 
Figure  16  and  the  symbols  have  meaning  as  defined  in  the  appendix. 
Figure  17  shows  the  mathematical  expression  for  two  elliptically 
polarized  plane  waves,  A(t)  and  S  (t),  traveling  in  directions  kx 
and  k2,  respectively,  with  a  time  phase  difference  f3(t): 


A(t) 


B(t) 


Re  ["(a  'u  +  A  e^1 

L\  e  9i  $ 

Ref(  B  u  +  B 
L\  e  6s  t 


2  )  ej(«.t  -  •  ?)]  (1) 

5  )ei(»t-ii  -?+9<t>]  (2) 


As  shown  in  Figure  16,  consider  a  pair  of  antennas  parallel  to  the 
x-axis  and  located  at  (x  =  0,  y  =  a,  z  =  0)  and  (x  =  0,  y  =  -a,  z  =  0), 
respectively.  The  complex  voltages,  and  V^,  induced  in  the 
antennas  are; 
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12  Expected  Sky-Wave  Modes  for  Saigon  in  June  1967 
Showing  Attenuation  in  dB  with  Respect  to  Stronger 


December  1967  Saigon  105^  106.  3^E  Sunspot 


June  1967  Churchill  58.  7N  94. 2W  Sunspot  No.  90 
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FIGURE  14  Expected  Sky-Wave  Modes  for  Churchill  in 
December  1967  Showing  Attenuation  in  dB 
with  Respect  to  Stronger  Mode  (U) 


December  1967  Churchill  58. 7N  94. 2W  Sunspot 


v  ~ 

a 


A(6j.  )  e^1  +  B(62)  e' 


j(f(62.)  -  f(5x  )  +  6)  j*s 


Vb  ~  A(6X  )  e"'^1  +  B(62)  ej(f(6s)  "  f(&1  *  +  e"^2 

The  net  result  may  be  interpreted  as  follows:  The  amplitude  variation 
of  A(6X  )  is  strongly  dependent -on  the  time  phase  bx  between  the  two 
orthogonal  components  Aq  ,  A$  ,  of  the  vector  X.  A  similar  remark 
holds  for  B(62).  The  variations  occur  not  because  of  the  polarization 
of  the  wave  but  because  of  the  changes  in  the  polarization.  The  phase 
variation  of  either  Va  or  V^  is  also  strongly  sensitive  to  changes  in 
61  or  62 . 


Mathematical  Description:  A(t)  +  B(t) 


Where 


1 

1 

m 


A(t)  =  Re(A  u.  +  A  ej  (,)  u  )  ej  <“l  ‘  ’  r) 

e  -9  4  -0 

B{t)  =  Re[ (B  +  B  ej6*(t)u  ,  j  e(t)]  eJ<»‘  -  $2  •  ?> 
0  4  0 


Polarization: 

Direction: 

Phase  Difference: 


In  general,  elliptical 
,  k2 

P(t) 


FIGURE  17  Wave  Interference  (U) 


1 


Figure  18  shows  a  phasor  representation  of  the  antenna  voltages  Va 
and  Vb,  This  figure  illustrates  how  the  magnitudes  of  Va  and  Vb, 
as  well  as  their  relative  phase  difference  y,  depend  on  the  magnitudes 
and  phases  of  the  interfering  waves:  A(  6X  ),  B(62  ),  p‘(t). 

As  p  1  ( t)  varies  randomly,  the  probability  distribution  of  the  angle 
that  Va  makes  with  the  real  axis  is  symmetrical  with  respect  to  ijn  . 
Therefore,  this  angle,  averaged  over  many  random  values  of  p‘(t) 
is  <!fx  .  A  similar  argument  holds  for  V^  The  total  average  phase 
difference,  y,  between  V&  and  Vb  is  therefore  equal  to  2\|fx  . 


This  last  statement  may  be  expressed  mathematically  by  the  theorem 
that  follows.  The  theorem  holds. for  all  the  distributions  (3'(t)  that 
allow  an  integration  with  respect  to  time  to  be  replaced  by  an  inte¬ 
gration  with  respect  to  (3. 

Theorem: 

Let 

Ej  =  Ae’*1  +  Be-'**3  +  =  |  E;  |  e^Yl 

=  Ae'ifc  +  Be'1**3  -6)  =  |  E»  |  eiv3 

Ys  =  2 


then 

Yi  " 


If  A  and  B  are  time  varying  but  |  A |  <  (  B|  at  every  instant  of  time, 
then  the  branch  point  never  moves  inside  the  contour;  hence  under 
these  conditions  the  theorem  still  holds.  From  this  we  may  con¬ 
clude  that  the  following  corollaries  are  true.  ■ 


Corollary 


If  A  and  B  are  time  varying  but  \  A(t)|  >  |  B(t)|  con¬ 
sistently,  y  -♦  2i|x  as  the  sampling  time  increases. 

s'  ~\ 


Corollary 


/ 


If  A>  B  some  of  the  time 
and  B>  A  the  rest  of  the  time 

but  |  A  |  >  |  B  j  ,  the  average  phase  y  lies  between 
2  0X  and  *+  \Jr3  ). 

Verified  by  computer  experiments  and  proven  analytically. 


PROBABILITY  DISTRIBUTION  OF  THE  PHASE  DIFFERENCE 
AT  THE  ANTENNA  TERMINALS  (U) 

Computer  analysis  furnished  information  as  to  the  probability 
distribution  of  y  as  p(t)  was  allowed  to  vary  randomly.  Two  main 
cases  are  of  interest,  but  for  all  cases  g(t)  assumes  an  independently 
random  variation. 


Case  I:  A  and  B  maintained  constant 

A  =  10,  B  =  1,  3,  5,  8,  10 

Case  II:  A  and  B  vary  randomly  with  Limits 

0  A  i  10;  0  s  B  <:  1,  3,  5,  8,  10 

Case  IIA:  Same  as  Case  II;  however,  two  distributions 
are  obtained: 

a.  A  s  B 

b.  B  s  A 

In  the  first  case,  2,  000  random  vaLues  of  gyieLded  a  probability 
distribution  simiLar  to  that  shown  in  Figure  19  for  alL  vaLues  of  B. 
Figure  20  represents  the  case  A  =  10,  B  =  5.  Two  conclusions 
can  be  drawn  from  this  portion  of  the  study.  For  each  distribution 
the  average  is  2^,  that  phase  associated  with  the  stronger  ray.  It 
is  aLso  noted  that  the  most  probable  phase  does  not  correspond  to 
either  of  the  two  directions  of  arrival.  Thus  under  the  condition  of 
two  mode  propagation,  with  stable  amplitudes,  the  phase  indicated 
most  often  is  in  error  with  respect  to  both  rays. 

The  second  case  considered  3,000  random  independent  samples  of 
A,  B,  and  gt  The  probability  distribution  shown  in  Figure  21 
illustrates  0  £  A  £  10;  0  g  B  s  5.  There  are  two  local  peaks  in 
the  curve  as  well  as  a  main  peak.  The  main  peak  corresponds  to 
2^,  while  the  other  two  correspond,  one  to  the  average  phase  of 
the  distribution,  the  other  to  2f-8,  the  phase  associated  with  the 
weaker  ray,  B.  It  is  noted  that  the  average  phase  no  Longer  cor¬ 
responds  to  2  ,  but  is  a  function  of  the  upper  limits  set  on  B. 

This  error  is  shown  in  Figure  22.  The  error  is  zero  in  the  limiting 
case  of  a  vanishing  weaker  ray  and  reaches  a  maximum  of 

Iq-L-la. 

2 

when  the  upper  limits  of  A  and  B  become  equal. 

The  final  case  is  but  an  elite  ration  of  Case  II.  For  each  calculation 
of  y»  the  computer  distinguished  whether  A  was  greater  than  B 
and  vice  versa.  In  this  way,  two  probability  distributions  were 
constructed,  the  sum  of  which  is  represented  by  the  previous  case. 
Figure  21  shows  the  distribution  in  red  for  B  >  A  and  in  black  for 


02J»V 


■i 
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A  >  B.  Notice  that  the  histograms  do  not  overlap  except  for y  =  19 
where  rounding  off  of  the  phases  to  the  nearest  degree  caused  a 
slightambiguity".  For  these  two  distributions,  the  averages  are 
2^  and  2ijfs,  respectively.  From  these  results  one  may  reach 
two  major  conclusions: 

(a)  When  one- ray  is  stronger  consistently  but  amplitudes  are 
varying,  averaging  still  yields  the  phase  corresponding 
to  the  stronger  ray. 

(b)  The  second  conclusion  is  that  if  one  could  differentiate 
as  to  which  ray  was  stronger  for  a  given  example,  then 
averaging  over  many  samples  could  be  performed  to  give 
the  directions  of  arrival  of  each  mode. 

Further  investigation  of  this  possibility  is  necessary. 

5.  EXPERIMENTAL  RESULTS  (U) 

Basically,  the  system  used  for  direction  finding  is  a  phase 
measuring  device.  To  apply  the  time  averaging  theorem  discussed 
in  the  previous  chapter,  the  phase  between  the  voltages  induced 
in  a  pair  of  horizontal  dipoles  must  be  measured  at  successive 
intervals  of  time.  The  heart  of  the  system  is  a  90-degree  hybrid 
which  is  conveniently  used  to  convert  the  problem  of  measuring  phases 
into  one  of  measuring  voltages. 


Vr 


Input 


V3  a 


90°  Hybrid 

Output 

V3 


V4 


FIGURE  23  Voltage  References  for  the  90-Degree 
Hybrid  (U) 
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This  particular  application  of  the  hybrid  was  developed  by  Jules 
A.  Cummins  of  the  Syracuse  University  Research  Corporation, 
and  is  given  here.  The  voltage  references  are  shown  in  Figure  23 


Let: 

Yx  =  Kl  =  |V1  | 


(3) 


V3  =  |'V3|  eJY  (4) 

Then  if  the  hybrid  is  ideal,  the  voltages  V3  and  V4  are  related  to 
Vx  and  V2  by: 


V3 

v4 

where 


-  ua 

41. 


(i-jqejY) 
(-j  +  qe^Y) 


q 


(5) 

(6) 


from  which  one  obtains: 


Kl8  - 

i  1  Vx  |  2  <l  +  q2  +  2q  sin  y) 

(7) 

K|3  = 

||V1  |  2  (!  +  q3  -  2q  sin  y) 

(8) 

Adding: 

Kl3  +  |V4)3  =  |  j  2  (1+q3)  =  IVJ2  +|V2 

(9) 

Subtracting: 


Kl* 

-  Kl2  = 

2  lVi  1  1  vs  1  sinY 

(10) 

|v,i*  -1 

Y*\ 

s 

(ID 

sin  y 

2  Vx 

v2 
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Subtracting  2  |  Vx  \  J  V2  |  from  both  sides  of  Equation  9: 


|v3|2+|v4|  =  -  2|V,|  |va|  =  (|vj  -  |V3|)2 

2lvll  |v,|  =  iv3|a  +|V4|=  -(Iv,!  -  |v2|)3 

(12) 

Now  Equation  11  can  be  written  as: 


sin  y 


(13) 


The  determination  of  y  involves,  therefore,  three  voltage  measure¬ 
ments.  The  usual  restriction  of  having  \  Vx  |  =  |V2|  for  this  type 
of  phase  meter  can  be  thus  removed  by  measuring  the  difference 
between  the  input  voltages. 

> 

The  entire  data  collection  system  is  shown  in  Figure  24.  To  simplify 
matters,  only  the  salient  features  will  be  discussed.  The  antenna 
consists  of  two  pairs  of  dipoles  in  quadrature,  with  apertures  of 
20  feet.  Either  pair  may  be  utilized  for  phase  measurements.  They 
are  matched  to  a  twin-channel  receiver  by  180-degree  hybrids. 


The  receiver,  a  Racal  P.153A,  is  phase  and  amplitude  matched  by 
manual  compensation  done  with  a  variable  delay  line  and  a  variable 
attenuator.  The  two  signals  from  the  receiver  are  branched  into 
two  paths  to  determine  the  variables  necessary  for  solution  of 
equation  13; 

a.  The  signals  are  detected  and  differentially  compared 

.  to  determine  |  |  -  |  V  2  |  . 

b.  The  signals  are  fed  into  the  90-degree  hybrid,  with  \  V3  | 
and  |  V4  |  being  detected  at  the  output. 

The  three  voltages  are  then  processed  by  analog-to-digital  converters 
for  recording  on  magnetic  tape. 
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The  serializer  in  Figure  24  locates  and  serializes  the  data  flow 
onto  the  magnetic  tape.  After  processing  data  fed  to  it  by  the 
analog- to- digital  converters  and  channeling  it  to  a  magnetic  tape 
-recorder,  the  serializer  interrogates  the  converters  for  new  data 
and  repeats  the  process  to  yield  a  new  phase  measurement.  Each 
cycle  requires  100  ms,  limited  by  the  time  constant  of  the  converters. 
Thus  the  system  measures  the  phase,  y»  at  a  rate  of  10  'samples 
per  second.  The  magnetic  tape  is  then  processed  by  a  Scientific  Data 
Systems  Model  930  computer  utilizing  a  computer  program. 

There  is  also  provision  on  the  serializer  for  manual  data  to  be  read 
onto  the  tape  prior  to  the  phase  measuring  process.  Manual  data 
entered  included  frequency,  time,  date,  and  ionosphere  height 
measured  by  an  on-site  ionosonde. 

Standard  procedure  was  to  measure  voltages  for  approximately  one 
to  two  minutes  and  thereby  enter  several  hundred  successive 
measurements  onto  tape.  The  computer  calculated  the  phases, 
printed  them  in  time  sequence,  and  averaged  them.  By  rounding 
these  phases  to  the  nearest  degree,  a  histogram  was  compiled  and 
printed  for  comparison  with  the  theoretical  distributions  presented 
in  the  last  chapter. 

This  process  was  done  twice  for  each  experiment,  first  in  the 
East-West  dipoles,  then  on  the  North-South  pair.  Thus  two  average 
phases  were  generated,  2 »!/x  and  2 tyy.  respectively.  Although  it  has 
been  assumed  in  the  previous  chapter  that  the  azimuthal  angle,  0, 
was  known,  the  two  average  phases  provide  both  0  and  0  of  the 
stronger  wave.  Referencing  0  to  degrees  east  of  north: 

2  TT  2L 

2 lit  =  — —  sin  0X  sin  &  (East-West  Pair) 

X  \ 

2  TT2t 

Z\)  =  —  sin  0X  cos  0X  (North-South  Pair) 

y  a 

Dividing  the  first  equation  by  the  second  and  solving  for  the 
azimuthal  angle: 
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The  elevation  angle  can  be  found  as: 


(2<sg2  +  (2*y)2 


(sin201 


+  cos2  0j_  ) 


0! 


",  1 


sin 


-i ;  - 


<2*/ 


+  <2'V>  J 


2na/\ 


:) 


(15) 


Experimentation  was  limited  to  CHU,  a  station  which  transmits 
time  signals  from  Ottawa,  Ontario,  Canada.  Located  nearly  due 
north  of  Syracuse,  CHU  transmits  a  vertically  polarized  signal  at 
7.  333  MHz  with  a  radiated  power  output  of  10  kilowatts.  The  re¬ 
sults  from  experimentation  with  CHU  transmission  are  shown  in 
Table  3,  where  they  are  compared  with  expected  results.  Since 
the  primary  goal  is  elevation  angle,  resolution  azimuthal  results 
are  not  included.  The  angles  of  arrival  of  the  IF  and  2F  modes 
are  expected  to  be  near  27  degrees  and  14  degrees,  respectively, 
assuming  an  F-layer  height  of  250  kilometers.  Also,  the  2F  mode, 
is  expected  to  be,  on  the  average,  18  dB  below  the  IF  mode,  cal¬ 
culated  using  National  Bureau  of  Standards  attenuation  procedures. 
Hence,  under  normal  conditions,  it  is  expected  that  averaging 
would  yield  the  stronger  ray. 


Arbitrarily,  relative  errors  in  angle  of  arrival  less  than  five  per 
cent  will  be  considered  satisfactory.  In  6l  per  cent  of  the  experi¬ 
ments  presented,  the  angles  measured  fell  below  this  limit.  The 
histograms  for  these  cases  were  very  similar  and  a  representative 
pair  from  the  March  27th  experiment  at  953  is  presented  in  Figures 
25  and  26.  The  phase  distribution  in  Figure  25  w as  obtained  from 
the  east-west  pair  of  dipoles.  Since  CHU  is  nearly  due  north, 
this  pair  is  most  sensitive  to  changes  in  d>.  The  average  of  this 
distribution  is  4.  381  degrees,  near  which  there  is  a  peak.  Also, 
there  is  a  peak  on  either  side  of  this  one. 


According  to  the  theoretical  computer  study  conducted  at  SURC, 
this  distribution  would  represent  amplitude  variations  with  the 
two  predominant  rays,  IF  and  2F,  arriving  via  different  azimuthal 
angles. 
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FIGURE  26  Experimental  Phase  Distribution  North-South  A 63 81 

Oriented  Dipoles  (U) 


Figure  26  shows  the  histogram  obtained  from  readings  on  the  north- 
south  dipoles,  which  are  phase  dependent  more  on  elevation  angles 
of  the  incoming  rays.  This  distribution  bears  ?.  marked  similarity 
to  the  theoretical  distributions  run  on  the  computer  lending  some 
validity  to  the  model  adopted  for  variations  in  h  and  3.  The  major 
difference  appears  to  be  the  more  Gaussian  appeara'hce  of  the 
experh-aenfs.!  hi.v'royvam.  A  plausible  expls ration  for  this  may  be 
found  in  the  theory  of  random  variable t.,  A 'v-.di.ng  lc  the  Ce  ntral 
Limit  Theorem,  the  more  random  variables  that  are  allowed  to 
affect  a  set  of  measurements,  the  more  Gaussian  will  be  the 
distribution  of  these  measurements.  In  this  problem,  although  only 
two  interfering  rays  are  considered,  many  more  are  actually  present 
at  the  antenna  much  weaker  than  the  two  primary  rays.  Their 
presence  would  cause  a  minimal  effect  on  averaging  for  the  stronger 
ray,  but  their  total  influence  would  be  to  produce  a  more  Gaussian 
appearing  phase  distribution. 

More  important  features  of  the  north- south  distribution  include  the 
pronounced  peak  near  the  average  of  21.  88  degrees  and  the  evidence 
of  a  weaker  mode  tending  to  widen  the  distribution  below  this  value. 
To  appreciate  the  value  of  time  averaging,  the  wide  variation  in 
measured  phases  should  be  noted.  An  instantaneous  measurement 
would  yield  an  uncertain  angle  of  arrival,  possibly  far  from  the 
angle  corresponding  to  the  stronger  ray. 

The  distributions  for  experiments  in  which  the  error  was  greater 
were  markedly  different.  The  test  on  22  March  at  1937  hours  is 
representative  and  the  histograms  are  given  in  Figures  27  and  28. 

A  general  characteristic  of  the  distributions  in  these  cases  is  their 
higher  dispersion,  which  indicates  more  severe  interference.  This 
is  especially  pronounced  in  Figure  28,  the  distribution  most  indi¬ 
cative  of  variations  in  elevation  angle. 

It  was  mentioned  previously  that  on  the  average,  the  2F  mode  would 
be  18  dB  below  the  IF,  thus  allowing  very  little  chance  for  h  to  be 
greater  than  one  under  normal  conditions.  However,  in  the  presence 
of  localized  ionospheric  absorption,  it  is  possible  that  the  IF  mode 
may  undergo  additional  attenuation  causing  the  2F  mode  to  interfere 
more  severely.  In  this  way,  many  values  of  y  may  be  computed  for 
which  the  2F  mode  is  stronger.  This  led  to  an  error  in  the  angle 
of  arrival  with  respect  to  the  normally  predominant  ray. 


52 


FIGURE  27  Experimental  Phase  Distribution  East-West  A6382 

Oriented  Dipoles  (U) 


Referring  to  Figure  28,  it  appears  that  the  2F  mode,  near  the 
peak  at  13  degrees,  figured  significantly  in  the  averaging,  which 
yielded  a  phase  below  that  associated  with  the  IF  mode.  Further 
evidence  is  afforded  by  the  many  values  of  phase  located  near  the 
peak  at  13  degrees. 

Some  experiments  indicated  an  angle  of  arrival  larger  than  that  of 
the  IF  mode.  Examination  of  the  ionograms  for  these  cases  often 
indicated  the  presence  of  an  Fx  layer,  which  does  not  normally 
support  a  reflection.  A  wave  reflected  from  this  layer  would  arrive 
at  an  angle  near  35  degrees  and  interfere  severely.  It  is  felt  that 
this  mode  was  present  in  these  cases  since  an  ionosonde  return 
indicates  probable  oblique  reflection  from  CHU. 

On  the  whole,  the  initial  experimental  results  are  encouraging  but 
much  more  investigation  is  necessary.  The  averaging  process 
yields  very  good  results  when  there  is  evidence  from  the  histograms 
that  the  IF  mode  is  substantially  stronger  throughout  the  sampling 
time.  In  these  cases,  there  is  only  a  slight  error  due  to  averaging 
values  of  y  for  which  the  2F  mode  becomes  predominant.  However, 
one  of  the  weaknesses  of  the  procedure  is  a  lack  of  reliability;  that 
is,  knowing  when  the  conditions  necessary  for  proper  application  of 
the  time  averaging  theorem  have  been  satisfied.  This  is  due  to  an 
inability  to  recognize  which  is  the  primary  ray  when  the  phase  is 
sampled.  This  will  be  studied  in  the  follow-on  work. 

TWIN  FRONT  END  ALTERNATIVES  STUDY  (U) 

The  purpose  of  this  study  is  to  investigate  the  tradeoffs  in  twin 
front  end  design  preference  based  on  stated  questionable  factors 
for  application  in  the  Amphibious  HF  D/F  System  being  developed 
at  SURC. 

The  front  end  is  to  have  simultaneously  tuned  channels  for  ampli¬ 
fication  of  directional  signal  components  in  the  high  frequency  range 
(0.  5  to  20.  0  MHz)  without  changing  their  relative  amplitudes  or 
phase. 

There  are  three  twin-front  end  system  alternative  approaches  as  shown 
in  Figure  29. 


Distributed 


Matched  Unmatched  Lumped  Front 

Channels  Channels  with  End  with  Matched 

Calibration  '  Channels 


FIGURE  29  Twin  Front  End  Alternatives  (U) 


a.  Twin  Front  End  Design  Preference  (based  only  on  operational 
desirability)  versus  Questionable  Factors  (U) 


The  following  table  states  the  twin  front  end  design  preference 
and  lists  the  questionable  factors  under  each  preference. 

1st  Choice  --  Distributed  Twin  Front  End  with  Matched  Channels 
Actual  operational  advantages 
Technical  feasibility 
Cost  and  time  constraints 
Manufacturers  acceptance  of  task 

2nd  Choice  --  Distributed  Twin  Front  End  with  Calibrated  Unmatched  Channels 

Actual  operational  advantages 
Technical  feasibility 
Cost  and  time  constraints 

3rd  Choice  --  Lumped  Twin  Front  End  with  Matched  Channels 
Actual  operational  advantages 
Cost  and  time  constraints 

b.  Evaluation  of  Twin  Front  Er.d  Techniques  (U) 

The  following  factors  will  be  considered  and  reported  on  as  to 
technique  and  technical  specification  for  each  of  the  three  system 
alternatives: 

(1)  Channel  amplitude  and  time  delay  tracking 

(2)  Sensitivity 

(3)  Dynamic  range 

(4)  Frequency  control 

(5)  Detection  and  signal  processing 

(6)  Spurious  responses 

(7)  Remote  programming 

(8)  Calibration 

(9)  Human  engineering 
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Special  attention  will  be  paid  in  reporting  on  the  following 
main  areas  of  interest: 

(1)  RF  Processing 

(a)  Input  Stages 

Shall  comment  on  devices  used. 

Possible  device: 

Parametric  up  converter 
Bi-polar  transistor 
FET 

(b)  RF  Input  Impedance 

Shall  comment  on  input  impedance  specification  and 
antenna  impedance. 

(c)  Noise  Figure 

Shall  comment  on  receiver  noise  figure  specification 
for  the  three  alternative  systems. 

(d)  Signal  Handling  Capability 

Shall  comment  on  signal  handling  depending  upon  in¬ 
put  devices  used. 

Shall  comment  on  weak  signal  optimization  versus 
strong  signal  tradeoffs. 

(e)  Tuning 

Shall  study  tuning  technique  to  be  used  and  time  re¬ 
quired  to  tune  the  receiver.  Possible  tuning  techniques 

Electronic  (voltage  tuned  variactors) 

Broadband 

Mechanical 

(f)  Spurious  Response  Suppression 

Shall  comment  on  technical  specifications  for  the 
following: 

Image  Rejection 

Intermodulation  products  (state  conditions) 

IF  rejection 
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(2)  Frequency  Control 

All  system  alternatives  will  be  capable  of  handling  a 
minimum  of  two  independent  signal  frequencies  at  the 
same  instant  of  time.  The  technique  for  frequency  control 
for  each  of  the  three  system  alternatives  will  be  studied. 

(a)  Digital  Synthesizer 
Phase  lock  loop 
Integrated  circuit  reliability 

10,  100,  1,000,  10,  000  Hz  channel  spacing 
Single  crystal  reference 
Remote  digital  control  available 
Reduced  spurious  responses 

(b)  Direct  Synthesizer 

High  speed  frequency  switching 
Interpolation  in  any  decade  increment 
Remote  control  available 

(c)  Hybrid  Synthesizer 
Low  power  consumption 
Maximum  flexibility 
Electro  nic  tuning 

Frequency  synthesizer  constraints  will  be  studied. 
(Injection  level,  impedance,  spectral  purity  for 
local  oscillator. ) 

Frequency  read-out  techniques  will  be  studied. 

(3)  Signal  Processing  and  Detection 

The  Twin  Front  End  system  must  be  capable  of  detecting 
all  types  of  modulation  inherent  in  the  HF  frequency  range 
of  0.  5  to  30.  0  MHz  including  AM,  FM,  USB,  and  LSB. 

The  last  intermediate  frequency  shall  be  no  lower  than 
455  kHz.  The  last  IF  bandwidths  should  be  800,  4,  000 
and  10,  000  Hz.  The  desired  shape  factor  of  the  IF  filters 
is  4:1.  (60  dB  BW:  6  dB  B\V) 


Linear  phase  characteristics  in  the  last  IF  are  desired. 
Selectivity 

Shall  comment  on  the  type  of  IF  filters  that  will  be  used 
and  the  technical  specifications  that  can  be  expected. 

Heterodyne  Scheme 

Shall  state  the  number  of  frequency  conversions  to  be 
used. 

Detectors 

Shall  be  capable  of  demodulating  the  above  mentioned 
modulations. 

AGC 

Shall  comment  on  the  following  parameters,  giving 
technical  specifications  where  possible. 

Dynamic  range 
Attack  time 
Decay  times 

AGC  technique  (IF,  RF,  LO  control) 

Miscellaneous 

The.  portion  of  the  Twin  Front  End  system  that  is  in  the 
shelter  will  meet  the  temperature  requirements  of  Mil 
E- 16400,  Class  IV. 

The  portion  of  the  Twin  Front  End  system  operated  out¬ 
side  of  the  shelter  will  meet  the  temperature  requirements 
of  Mil  E- 16400,  Class  I. 

Packaging 

Shall  comment  on  the  packaging  techniques  that  will 
be  used  in  system  fabrication. 


Shall  house  the  first  stages  of  the  front  end  at  the 
antenna  array  in  the  distributed  Twin  Front  End 
case.  The, last  IF  stage  will  be  the  first  subsystem 
stage  in  the  shelter. 

Summary 

Shall  report  the  above  parameters  and  techniques 
for  each  of  the  three  system  alternative  approaches. 
Specifications  and  block  diagrams  of  suggested 
systems  will  also  be  furnished  in  the  study  report. 
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Section  III 

FOLLOW-ON  WORK  PLAN  (U) 


Demonstration  of  feasibility  in  the  three  risk  areas  now  permits  a  rela¬ 
tively  straightforward  approach  to  the  accomplishment  of  formulating 
quantitative  specifications  for  a  complete  improvement  of  the  experimental 
system.  The  experimental  system,  as  now  envisioned,  will  include  an 
array  capable  of  narrow  beam  D/F  on  ground  waves  and/or  single  site 
fix  on  complex  sky-waves.  This  array  will  consist  of  optimal  feed  net¬ 
works,  preamplifiers,  and  a  Twin  Front  End  with  precision  matching. 

The  last  section  of  the  array  will  consist  of  electronic  processing,  analog- 
to-digital  converters,  and  an  RF  calibration  capability. 

The  shelter  will  accept  the  digital  data  inputs  from  the  array  and,  utilizing 
near -real- time  computer  analysis,  will  retrieve  fix  and/or  D/F  infor¬ 
mation  automatically.  This  information  will  be  made  available  to  the 
operator  through  computer  print-out  and  visual  displays.  Other  infor¬ 
mation  generated  includes  polarization,  frequency,  demodulated  signal, 
modulation  mode,  site  assessment,  D/F  fix  classification,  and  ionosonde 
information. 

This  follow-on  work  covers  completion  of  the  array  portion  of  the  system 
just  discussed.  Specific  tasks  to  be  accomplished  under  this  effort  are 
as  follows: 

a.  An  experimental  program,  involving  testing  of  active 
elements  with  optimized  feed  networks  and  preamplifiers, 
will  be  conducted  in  order  to  determine  the  final  array 
configuration.  This  program  will  include  computer  analysis 
and  propagation  testing. 

b.  In  order  to  determine  optimal  time  convergence  of  sampling, 
sampling  rate,  time  per  unit  sample,  and  fade  spectrum, 
further  propagation  tests  will  be  run.  These  tests  will  also 
add  additional  range  and  seasonal  propagation  information 

to  the  propagation  data  base  already  obtained. 

c„  Based  on  the  information  obtained  in  the  preceding  task,  array 
processing  electronics  will  be  designed  and  fabricated. 

d.  A  developmental  test  model  of  the  twin  front  end  electronics 
will  be  designed  and  fabricated  for  use  in  the  array. 
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GLOSSARY  (U) 


1,  2 


Re 

Im 


V  Be 

-» 

ki,3 

k 


COS  (X\Z 


cos  0  1>s  ,  cos  y1>3 

U„,  “al  U  ,  U  ,  U 

—0  -~<b  —x  —  y  — z 


Jl,3 


3(t) 


r=xu  +  yu  +zu 
—x  — y  — z 


Subscripts  1,  2  refer  to  waves  A  and  B, 
respectively. 

The  real  part  of. 

The  imaginary  part  of. 

The  projection  of  the  magnitude  of  the 
electric ormagnetic  field  vector  associated 
with  wave  A  propagating  in  direction  kj,. 

Corresponding  definition  for  wave  B. 

k(COSai.2  U  +  COS  01  3  U  +  COS  V-i  p  u  ). 
’  —x  »  — y  —z 

2tt  l\,  X  being  the  free  space  wavelength. 

The  direction  cosine  with  respect  to  the 
x-axis. 

Corresponding  definition  for  the  other  axis. 

Unit  vectors  directed  toward  positive 
direction  of  the  subscripted  coordinate. 

The  time  phase  difference  between  A^ 

and  A,  or  B„  and  B^.  1 

<b  0  d> 

The  time  phase  difference  between  the  two 
waves. 

The  radius  vector. 


*1,3  =kasin0i>3  cos§i^2 
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